Equilibrium-like fluctuations in some boundary-driven open diffusive systems 
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There exist some boundary-driven open systems with diffusive dynamics whose particle current fluctuations 
exhibit universal features that belong to the Edwards-Wilkinson universality class. We achieve this result by 
establishing a mapping, for the system's fluctuations, to an equivalent open -yet equilibrium- diffusive sys- 
tem. We discuss the possibility of observing dynamic phase transitions using the particle current as a control 
parameter. 
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I. INTRODUCTION 

In this work we consider systems of interacting particles undergoing diffusive dynamics, such as the Symmetric Simple 
Exclusion Process (SSEP) |H |^, |H 01 . Such systems, that can be described by the theory of fluctuating hydrodynamics, a 
coarse-grained description in terms of continuous degrees of freedom living in a continuous space, have already been the 
subject of intense investigation. For instance, Bertini et al. JH @, 0, [ft 0] have relied on fluctuating hydrodynamics to provide a 
quantitative analysis of large deviation properties of diffusive systems taken out of equilibrium by means of a boundary drive. 
Among the various large deviation properties investigated so far, those of the particle current play a special role. Indeed, current 
fluctuations have been known for a long time to be a central quantity since the work of Einstein II 1 OL ll ill which established that 
in equilibrium the current variance is proportional to the diffusion constant; current fluctuations characterize the likeliness of, 
and quantify, the system's excursions out of equilibrium. In the last decade, generic properties of these large deviation functions 
were discovered such as the fluctuation theorem which determines how the lar ge d eviation function of the current is changed 
under time reversal 

02BQ[H[Ii[I3[H[I!|23]. Parallel a PP roaches JUTil HI Hi HI Si which have been employed 
have revealed the possibility of new types of phase transitions where the particle current plays the role of a control parameter. It 
must also be mentioned that some results regarding the particle current statistics originate from exact solutions; this is the case 
for the totally asymmetric exclusion process Il27tl -a version of the SSEP where the particles' motion is strongly biased- or for 
the SSEP 112811 . with periodic boundary conditions. 

Fluctuating hydrodynamics not only encompasses interacting particle systems, but also applies to models involving at 
the microscopic level already continuous degrees of freedom, such as the Kipnis-Marchioro-Presutti l29ll (KMP) model of 
interacting harmonic oscillators that has served as a testbench to investigate the statistical properties of heat conduction. Both 
the SSEP and the KMP models will be the main focus of our efforts in the sequel, though we shall strive to keep our discussion 
general when possible. 

Our central motivation is to investigate the role of finite-size effects in one-dimensional open driven diffusive systems. We 
have been inspired by the results of Appert et al. ll28ll and by those of Derrida and Lebowitz ll27ll . In the former universal 
properties were seen to emerge in the statistics of the particle current for an equilibrium diffusive system with periodic boundary 
conditions, with the possibility, for certain classes of diffusive systems, to exhibit a current-driven dynamic phase transition. In 
the latter, where mutually excluding particles are subjected to a bulk electric field that drives the system far from equilibrium, 
universal features belonging to a different universality class have also been observed. Besides, Bodineau and Derrida J25l l26ll 
have shown that for a weakly bulk and boundary driven SSEP, a dynamic phase transition takes place. 

In the present work our interest goes to open systems, maintained out of equilibrium by putting them in contact with particle 
reservoirs at unequal chemical potentials, but the bulk dynamics itself remains reversible. Thus the nonequilibrium nature of 
our systems does not arise from an external bulk field but only from a boundary drive. We ask the following questions: (i) Do 
universal features in the particle current appear in an open system? If so, do they depend on the system being possibly driven 
out of equilibrium by a chemical potential gradient? (ii) Is the existence of a particle current capable of inducing a dynamic 
phase transition? 

Before entering the technicalities of our work, we would like to phrase the answers we have come up to issues (i) and (ii). To 
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question (i) we have the partial answer that at least for a class of systems -to which the SSEP and the KMP model belong- the 
current distribution does indeed display universal features. The latter do not depend on the system being in or out of equilibrium, 
and, quite remarkably, they are the same for an open system as for a closed system rf28ll . They belong to the Edwards-Wilkinson 
universality class 13011 . To question (ii) the answer is not straightforward: we have found at least a family of physical systems 
(among which the SSEP) in which the current large deviation function displays some singularity, indicating the existence of a 
dynamic phase transition depending on the scaling of the current one forces through. 

We shall begin in section [TT] by recalling what is known on the statistics of the current in an open boundary driven diffusive 
system. Section [III] is devoted to careful analysis of finite-size effects leading to establishing that in some cases fluctuations 
exhibit universal features. This section is supplemented by an appendix that describes the cases of the SSEP and the KMP model 
in detail. Our conclusions and yet open problems are gathered in section HVl 



H. CURRENT LARGE DEVIATIONS IN DIFFUSIVE IN BOUNDARY-DRIVEN OPEN SYSTEMS 

We consider a one-dimensional lattice with L sites, whose state at time t is characterized by the local numbers rij (t)'s (which 
may be discrete or continuous variables), j = 1, . . . , L. Our starting point is the assumption that, in the large t and L limit, with 
t/L 2 fixed, there exists a Langevin equation for a density field p(x, r) = rij(t'), with x = j/L and r = t 1 / L 2 defined over 
x G [0, 1] and t € [0, t/L 2 } which evolves according to 

d T p = d x (D(p(x, r))d x p{x, t) - £(x, r)) (1) 

where the Gaussian white noise £ has correlations (£(x, r)£(x', r')) = a ^?' T " S(x — x')S(t — r') that decay to zero as the 
inverse system size. The phenomenological coefficients D(p) (the diffusion constant) and a(p) depend on some of the details 
of the underlying microscopic dynamics. The system size L and the observation time t are large with respect to microscopic 
space and time scales. The system is in contact at both ends with reservoirs that fix the value of p at all times to be po at x = 
and pi at x = 1. 

Our interest lies in the statistics of the total particle current Q(t) accumulated up until time t and its large deviation properties, 
which, in terms of the field p, is formally expressed as 

r i f t/L 2 

Q{t) = L 2 Ax dT(-D{p(x,T))d x p(x,T)+Z(x,T)) (2) 

Jo Jo 



Our purpose is to determine 



= lim lnProb{Q(t) =J t} (3) 



or, equivalently, its Legendre transform xp(s) = max — sj} that can be obtained from the generating function of Q as 
follows 

m / e -sQ(t)\ 

V>(s) = lim — i '- (4) 

t— >oo t 

Using the Janssen-De Dominicis formalism [ f3lt] . we see that the generating function (e~ s< 2(*') can be cast in the form of a 
path-integral over two fields, 

where the action S is given by 

S[p,p]= I dx / dr pd T p + D{p)d x pd x p - ^-(d x p - X) 2 - \Dd x p 
Jo Jo 1 * 

where A = sL, and the path-integral runs over functions verifying the boundary conditions p(0, r) = po, p(l, r) = pi, and 
p(0,r) = p(l,r) = 0. It is very clear from the expression of the noise in (HJ or from the path-integral <(3j that a semi- 
classical-like expansion is valid in the weak-noise limit, which, translated in our language, is synonymous for a large system-size 
expansion. In short, the path-integral (0 is dominated by the saddle point of S. The reader is referred to Kurchan's lectures rf32"ll 



e -»ww) = ; VpDpz-^ p ^ (5) 



1 r t/L 2 

(6) 
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for a pedagogical account exploiting this language. We first change the response field into p{x, r) = p(x, r) — Xx. We write 
the saddle-point equations and we assume the saddle is reached for time-independent profiles p c (x), p c {x). Sufficient conditions 
under which this is so have been discussed by Bertini et al. |9]. Assuming this is indeed the case, the saddle point equations 



§ = 0,§ = Oread 



d T p = d x (Dd x p)-d x (ad x p), -d T p = d x {Dd x p) + y {d x pf (7) 
which, assuming a stationary solution, lead to 

D 2 ( Pc ){d xPc f = Kl + K 2 a(Pc), d xPc = D (PJ d fc + K i (8) 

where K\ and K 2 are A-dependent constants. With these equations, one may verify that the action evaluated at the saddle reads 
S[p c , p c ] = (t/L 2 )K 2 /2. We shall denote by p(X) = -Kijl (our definition of \i differs from that of 00 by a factor l/L). 
With these notations, for the large deviation function introduced in (0), we thus have ijj(s) — i^J^i. In practice, to determine 
p(X) explicitly, one must solve the differential equations ([8]l and fix the constants K\ and K 2 by means of the appropriate 
boundar y co nditions. A few comments are in order: these results are not new and they were first derived by Bodineau and 
Derrida 12411 . We propose as an illustration the explicit expression for /i(A) when the diffusion constant D is independent of the 
local density (we take D = 1) and when the noise strength cr(p) is a simple quadratic function. For a{p) = c 2 p 2 + C\p we find 
that (see appendix A) 

MM - I + l( arcsin y^) 2 foru> < W 
where u>(\, po, pi) is the auxiliary variable given by 

w(A, p ,px) = - e ClA / 2 ) U{ Pl - e- ClA /2 po) _ ^-^2 _ l)p \ (1Q) 

For the SSEP, a(p) = 2p(l — p) and one recovers the known J33ll34ll result (the notation z = e~ A is used in the formula (2.14) 
of JUt]), namely 

w(A, po, Pi) = (1 - e A )(e- A Po - Pl - ( e - A - l)p oPl ) (11) 

Another solvable model is the KMP chain of coupled harmonic oscillators, for which D = 1 and <j(/o) = 4p 2 (for KMP, p stands 
for the local potential energy field), for which we also have (0 but where the variable ui is now given by 

w(X, po, Pi) = ^(2(^0 - Pi) - 4A/30P1) (12) 

The case c 2 > and c 2 < are qualitatively different. In the latter, p(X) is defined over the whole real axis and is unbounded 
from above, while in the former p(X) is defined over a finite interval of A whose ends correspond to infinite currents produced 
by the build-up of infinite densities. For example, with c 2 = 4 and c\ = 0, that is for the KMP model, —5^7 < < 2p^- 

Finally, it is important to realize that p(X) is the leading order term in a large system-size series expansion. The origin of 
finite-size corrections is twofold. Of course there will be finite-size corrections arising from integrating out the modes describing 
fluctuations around the optimal profile {p c , p c }. However, fluctuating hydrodynamics, by definition, is unable to capture the 
details of the microscopic systems it describes. It must therefore be expected that model-dependent finite-size corrections will 
also emerge. We now proceed with determining the finite-size contribution of fluctuations around the saddle of the action (O 
within the framework of fluctuating hydrodynamics (that is, temporarily omitting contributions arising from the underlying 
discreteness of the lattice). 



III. FLUCTUATIONS AND UNIVERSAL BEHAVIOR 
A. Evaluating a determinant 

As in any saddle point calculation, we obtain finite size-corrections to the leading order result (e~ s ®) ~ e M ' si - ) * by in- 
troducing, in the path-integral (|5j, the fluctuations around the optimal profile p c and p c : 0(x, t) = p{x,r) — p c (x) and 
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(x, t) = p(x, t) — p~c(x). Then we expand the action © to quadratic order in and i 

[l(\)t f ,,_( 7a , r-Q 7a , , n/a - , a , n/a - a Ij. , D " . 



L 2 



J dxdr (^j)d T (j) + Dd x (j)d x (f> + D' d x p c 4>d x 4> + D' d x p c d x (jxj) + ^-d x pd x 



- a'd x p c <t>d x $ - G —{d x p c ftf 



(13) 



where D, a, and their derivatives with respect to the density, are evaluated at p c (x). The goal is to integrate out the quadratic 
action ( fT3l with respect to the fields <j> and <j>. This is the procedure that was followed in l28ll and that we carry out here as 
well. However, unlike the case of periodic boundary conditions dealt with in l28ll . in the present case, the quadratic action is not 
readily diagonalizable for its coefficients are space-dependent. It so happens that for one particular family of models, those for 
which D(p) is constant and a(p) is quadratic in p, this can actually be achieved. This remains a nontrivial task, given that the 
quadratic form to diagonalize in (foi l still possesses space-dependent coefficients. We have not been able deal with arbitrary D 
and a. 

B. Constant D and quadratic a 

We specialize the action $13[ to a constant D and a quadratic a. After performing the change of fields 

= (3e/5c)~V + d x p c $, 4> = d x p c i) (14) 
we note that ( fT3l , after tedious rearrangements, becomes 

M(A)i 



+ J dxdr \J>d r il> + Dd x i>d x ilj - fi(\)(d x rpy - —rj (15) 

It is remarkable that (TT3T > is now a quadratic form that can be diagonalized with standard stationary waves {suiqx} q with 
Fourier modes indexed by q — irn, n G N*. By comparison to (fT3b . we can interpret (T5[ as being the action corresponding to 
an equilibrium open system whose current fluctuations we study as a function of the conjugate variable Hifid Performing the 
change of variables (fT4l has allowed us to map the fluctuations onto those of an open system in contact with two reservoirs at 
equal densities. 

After integrating out the ip and ip fields, one arrives at 

VU*) = \p(sL) + ^-T f^-p( s L)) (16) 



' 8L 2 \2D 2 

where the FH index stands for "Fluctuating Hydrodynamics" where the function T has the expression 



Tin) = -4 

q—riTT, n> 1 



(qyjq 2 - 2u - q 2 + u) (17) 



Equation ( TToT l is the first new result of this work. It indicates that for systems whose fluctuating hydrodynamics description 
relies on a constant diffusion constant D and a quadratic noise strength a(p), current fluctuations involve a universal scaling 
function T. It is remarkable that exactly the same function T has appeared in the study of current fluctuations in closed systems 
in equilibrium, with a different scaling variable though. As can be seen from its explicit expression ( fTTb , the scaling function T 
has a singularity when its argument approaches 7r 2 /2 from the right real axis. In II28L 13511 this was interpreted as the presence of 
a first-order dynamic phase transition, for systems with periodic boundary conditions. In the present case, this also opens up the 
possibility of a similar phase transition on condition that there exists a regime of A for which 

tpW > ( 18 ) 



2D 2 ' 

Before we discuss whether a phase transition can indeed occur, we must address another pending issue. 

C. Microscopic details matter 



The expression for ip m {\/L) = j^p(\) + gj^J 7 ^^j5-m(A)J obtained from fluctuating hydrodynamics ignores the possibility 
that finite-size corrections of the same (D(L~ 2 ) order as the universal corrections will appear when one relies on the original 
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model defined on a lattice. In the appendix, which is based upon methods developed by Tailleur et al. I36ll37ll . we are able to 
evaluate the contribution of lattice effects for two specific models. We show that for the SSEP and for the KMP model they do 
introduce 0(L~ 2 ) terms that add up to the universal contribution found from fluctuating hydrodynamics. Let us look into those 
microscopic details more precisely, first for the SSEP, then for the KMP model. 

The open and driven SSEP consist of particles hopping to either of their nearest neighbor sites on a lattice of L sites, in contact 
with particle reservoirs connected to sites 1 and L. Particles are injected into site 1 (resp. L) with a rate a (resp. S) and are 
removed from site 1 (resp. L.) with rate 7 (resp. f3). These reservoirs impose densities po = and p\ — at sites 1 and 
L. While in the fluctuating hydrodynamic formulation the reservoirs enter current statistics through po and p\ only, when one 
wishes to capture phenomena beyond leading order, lattice effects and microscopic details start playing a role. For the SSEP, as 
presented in appendix B, introducing the auxiliary constants a = and b = -g^j, we find that 



</>(») =*.(») - £2 MA) + C(L- 3 ) 

=r"< A > - t"< a > + & (^"< A >) +otL ^ 

Note that this result is compatible with the exact expressions of the first three cumulants of the current obtained in ll33ll . 

The KMP model is also a lattice model in which L harmonic oscillators whose positions xj are coupled (we use the Ito 
convention and the Giardina et al. §5$\ version of the KMP model) 

dx 

2<i<iV-l, -j^- = -Xj +x j+ irjj >j+ i - Xj-irjj- ltj (20) 
and the chain is in contact at both ends with heat baths imposing temperatures T\ and T^, 

= - (jl + X X - + £2771,2, ^ = - Q + 7i) X L + ^2 lL T L ^ L - X L -lVL-l,L (21) 

where £1, ^, (for 1 < j < L — 1) are Gaussian white noises with variance unity, and 71, jl set the time-scale of the 

energy exchange with each reservoir. We refer the reader to Giardina et al. JUlEsl] for further details and connections between 
the SSEP and KMP. It is also shown in appendix B that for the KMP model we have 



i>{a) = j-p(X) Z71 ^ MA) + — T I ^2M(A)J + 0(L-*) (22) 

Both ( TT~9b and d22b reveal that taking into account microscopic details of the systems leads, as expected, to nonuniversal correc- 
tions to the current large deviation function. Whatever the form of these nonuniversal contributions to tp(s), it can be seen that 
the relevant piece of information regarding the possibility of a phase transition is contained in the universal part of ^ FH (s). 



D. Is a dynamic phase transition possible? 

For systems having a constant diffusion constant D (which we set to D = 1) and a quadratic a(p) = c^p 2 + c\p, the explicit 
expression of p,(X) obtained in © allows us to probe the criterion ( TT~8T > for the existence of a phase transition. Working at fixed 
A = sL. no first-order phase transition can occur because the condition (fT~8T>. or equivalentlv, u(X. po, pi) = 1, cannot be fulfilled 
on the real axis of A. In the original variable s, however, things are different, since in the large system size limit, and for C2 < 
only, the singularity in the complex plane of s eventually hits the real axis at s = 0. To be more explicit, using ([Toi l, one notices 
that for c 2 < 

A ^00, ^(A)~-^A 2 (23) 
2c 2 

so that, after inserting into (O and taking the asymptotics, one arrives at 

,. V'Q(s) 2c 1 2 . c l I ,3 , / 3n , 0/n 

lim — - — = s +— \s\ + o{s ) (24) 

L^oo L C2 OTT 

The singularity at s = reflects the existence of a dynamic transition in terms of the total particle current, but of higher order. 
The same transition existed for systems with periodic boundary conditions (see (62) of [,28]) and was noted earlier by Lebowitz 
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and Spohn (in (A. 12) of fl5tl ). The effects of this transition can be seen on the correlation functions {HI which become 
long-ranged. Note also that in this scaling limit (l24l i does not depend on the reservoir densities anymore, because the op- 
timal profile able to carry such large currents settles to density 4, but in vanishingly small region around the system's boundaries. 

For systems with c 2 > 0, that is systems with attractive interactions, such as the KMP model (for which a(p) = 4/? 2 ), no 
phase transition can be observed, but the trivial one occurring at infinite densities (akin to a Bose condensation). There exists a 
set of numerical simulations by Hurtado and Garrido i40ll for the KMP model which actually confirm that no phase transition is 
observed. This negative result is in contrast with -but does not contradict- that of Bertini et al. (|9l 12811 in which it was shown 
that a phase transition exists, for periodic boundary conditions, when a" > 0. 



IV. OUTLOOK 



We have shown that in a family of diffusive systems driven out of equilibrium by a chemical potential gradient, the total 
particle current exhibits universal fluctuations. These belong to the Edwards-Wilkinson universality class and they are of the 
same form as that previously found in closed equilibrium systems. Our results apply to diffusive systems characterized by a 
constant D and a quadratic a. We have used a mapping of the system's fluctuations to those of an equivalent open system in 
equilibrium. We have hints that this mapping can be extended beyond quadratic fluctuations: for the SSEP, we can actually 
prove that a similar mapping applies to the full process ll42ll . Our main concern lies in that our results are indeed limited to the 
case D constant and a quadratic. It would be of great interest to find out whether similar universal properties hold for generic 
D and a. Perhaps is this a fortuitous coincidence, but Tailleur et al. IL36L l37ll ran into similar restrictions when mapping the 
density profile large deviations in boundary-driven diffusive systems onto their equilibrium counterparts. We see here subjects 
for future research. 
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Appendix A 

In this appendix we prove (19110b . We assume that 

D{p) = l ; a{p) = c 2 p 2 + cip 



(25) 



with ci > and the boundary conditions p(0) = po, p(l) = p%. In order to find the explicit expression of /i(A) we start from 
the implicit equation found in ll24ll specialized to D = 1, namely 



MA) = -K 



P1 h 1 
ap—F= 



with K determined by 



We know that the optimal profile verifies 



rpi i 
A = / dp- 

PO ° 



- 1 



d x p = qVl + 2Ka 



the solution of which takes the form 



provided /, #o and 9\ verify 



p{x) = -- + /sinh [2(6 + (9 1 - 6 )x)] 



(26) 

(27) 

(28) 
(29) 



h - do? = ^c 2 Kq 2 



j2 = 2c 2 ~ C 1 K 



4c 2 2 K 



(30) 
(31) 
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and the boundary conditions p(0) = po, p(l) = p±. One performs the change of variable 



dx= 1 - - AP (32) 
q VI + 2AV 



in J2"6l > and (|271 l, This yields 

2 



M (A) = (0i-0 o r (33) 

C-2 



2 cip cosh26' 1 - ^Jc\ +Ac\f 2 p sinb^x 
A = H m - — (34) 

c i ci/9icosh20 o - V c i + 4c|/ 2 pi sinh2# 



where we have used OTb to eliminate K in favor of /, together with the boundary conditions 



Po = + f sinh 29 , Pl = -iL + / sinh 26 1 (35) 

2c 2 2c 2 



We are left with eliminating /, 9q and 6\ from (|33H35b . Isolating first f 2 from ( f34-b one obtains 

c 2 z 2 po - 2zp pi cosh 2(6*1 - #o) + Pi 
4 C 2 (z po sinh 26*i — pi sinh 26q 



* = ~ ; . ; ITT2 — (36) 



where we have set z = e ClA / 2 . Grouping / 2 with the square at the denominator in (l36l l. one eliminates / using the boundary 
conditions (f35T >. This enables us to isolate cosh2(0i — 9q) in d36t and one gets sinh 2 ((9i — 8q) = lu with 



w = -§ (l - z 1 )(ci(pi - zp ) - c 2 (z- l)poPi) (37) 

and finally 

r-|(axg S inhV7J) 2 2 foro,>0 (3g) 
^+^-(arcsin v^w) foiuj < 

In the limit c\ — > which is relevant for KMP, cj becomes 

w = ic 2 A(2(p - pi) - c 2 Xp Q pi) (39) 



Appendix B 

A. Simple Symmetric Exclusion Process 

We consider a SSEP on a one-dimensional lattice with L sites, in which particles are injected to the leftmost site j = 1 (resp. 
rightmost site j = L) with rate a (resp. 5), and removed with rate 7 (resp. /3). The master operator governing the evolution of a 
microscopic configuration of occupation numbers {%}j=i,... l can be written in the form 

i<fc<L-i (40) 
+ a(e~V+ - (1 - fn)) + 7 (eVf - ni) + <5(e"V+ - (1 - n L )) + /3 (eV^ - h L ) 

In ( f4Qb we are using a spin basis: the eigenvalue of the Pauli matrix cr| is 1 if site j is occupied, and —1 if it is empty (fij = 1+ ° 3 
has the eigenvalue rij = or 1). We now remark that 

e s ^=i ww(A)e- s ^=i^ = W L (s(L + 1)) (41) 
where Wl is the operator counting the total current across site L only, the expression of which reads 

Wi(s') = XI ( CT fe CT fe+i + CT fe °fc+i -n*(l-n fe+ i)-nfc + i(l-n fc )) 



K/c<L-l 



+ a(af - (1 - ni)) +7 (ai - ni) + <5(e s '<7+ - (1 - n L )) + /3 (e s 'tT L - n L ) 



(42) 
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with s' being conjugate to the time-integrated current through site i = L. Owing to fiTt . The largest eigenvalue of Wl(.s(L + 1)) 
is ^p(s), that is the largest eigenvalue of ( |40b . We use, for each lattice site, a Holstein-Primakoff like representation 0711 



cr+ = 1 - F + F+ - 2FF+ + F 2 F+ 1 a~ = F - F 2 F+ 
which also leads to h = F + FF + — F 2 F + . The bulk contribution to the evolution operator (|42| | now reads 

W L , bulk (A) = - ]T ((F j+1 FMF+ x - F+) + (F i+1 - FtfFfFf^) 



(43) 



(44) 



We represent e w ^ s ')* by means of a path-integral 0~7i l4lll involving coherent states related to the operators F and F + , which 
we shall denote by 0(t) and 0(r). This leads to an action 



= / di 



Si .bulk [? 

while the boundary terms are given by 



L L-l 

0^0? + - <t>j)(<i>j+i - 4>j) + (0j+i - <i>j) % <i>jh+i] 

3=1 3=1 



(45) 
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L, boundary [ 



di [a0i - (a + 7)<£i0i 



(46) 



dt 



-e- s + ((e- s - 1)0 L + 1)0 L + 1 



Since we are interested in the large-time behavior, we shall proceed with a saddle point approximation at fixed t but as 
L — > oo, keeping the system size L fixed (this is possible due to our saddle-point equations being stationary). We use the 
notation Vj0 = 0j+i — 0j. The saddle point equation obtained by differenting Sl with respect to 0j reads 



(Vj-0 + 2V i 00 j 0j+i) - (Vj-10 + 2V j _i00j_i0 i ) = 



and thus there exists if i such that 



-V J + 2gi 
V,0 = - 



203 

Writing the variational equation with respect to <frj and using (|48| ), we obtain 

|±l+^zl (4^-0 J 2 + J . + i0 3 ._i) = O 
03+10J03-1 

which we multiply by 

03+1 " 03-1 7 



y 3 + l fj-1 



so that 



4^i 2 - (V,- 



4if2 - (V 3 -_i< 



which leads to the existence of another constant K2 such that 

^2 a 1^2 



(V J -0 J -) 2 =4^ + iT ; 



2<Pj(Pj + l 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



We thus obtain that when evaluated at the saddle, S[<f>, 0] = — £-^— -K%. Besides, it is possible to solve the bulk saddle point 
equations d48l ) and d52l : 



2 < j < L - 1, 0j ; = -Asmh[(j - l)B + C], <j) 3 ,= E + — tanh ^ ^ 



(53) 
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where A, B and C are related to Ki and K2 by if 1 = — =j sinh-B, A'2 = 4sinh 2 -j. At this stage we write the saddle point 
equations corresponding to the fields located at the boundaries j = 1 and j = L. At j = 1 this leads to 

= 2 - 01 + 2(0 2 - 0i)020i - oc(f>i - 70i (54) 

and 

= 02 - 0i - (02 - 0i) 2 02 + a{l - 0i) - 70i (55) 

This immediately sets the constant E appearing in d53l ) to E — = po, and further imposes that a = = fj^-g- Due 
to the latter relation between B and C, only two unknowns A and £? remain to be determined. This is done by writing the two 
saddle point equations at j = L and by substituting the solution 1531 . The additional constraints on A and B (or C) are 

A[(z - l)p + 1] [zp (pi - 1) - P0P1 + pi] ' 

and 

sinh[(£ - 1)S + C + e] + - sinh C = (57) 

a 

where sinh 2 | = ui, p\ = -J^ is also the density at site L and where the variable lu is exactly that defined in (fTTb with s' instead 
of A. Finally, we eliminate C to obtain B as the solution to 

sinh 2 [{L - l)B + e} = (a 2 + b 2 + 2a6cosh [(L - 1)B + e] ) sinh 2 B (58) 

Equation (l58l can be solved in powers of 1/L: to leading order, B and C are 0(1/ L), while e is 0(1), and thus one has 
B = j^e = \ arcsinh y/tH. To the next order one has 

^l(s') =^-(-1 + Vl + a 2 sinh 2 B) + ^-(-1 + Vl + 6 2 sinh 2 B) 

+ (L - 1) sinh 2 ^ (59) 

^-^^MOO + 0(L-») 
This proves the result announced in ( fT9l . 



B. Kipnis-Marchioro-Presutti model 

For the KMP process, one writes a Langevin equation for £j = \x 2 based on (1201) . Using the Ito discretization scheme, this 
leads to 

-rr=3i- Ji+i (60) 
dt 



where the local energy current is = ei — £; + i + 2 v ^7e7+i?7j,i+i (1 < i < i — 2), and ji = 71 T\ — + 2V7i/7T£i> 

(61) 



Jl+i = —7lTl + %7l£l + ^Vtl^l^l- Using the Janssen-De Dominicis formalism, one is again led to 

( e - sQ ) = J v-;n-,c s ■ 

where the action has the expression 



J 3=1 J j=i 

+ 2 7l J dt [-Ti(ei - s)((ei - s)ei + 1/2) + - s) £l ] 
+ 2 7L J dt [-T L {e L + s){{e L + s)e L + 1/2) + (e L + s)e L ] 



(62) 
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With the change e'- = Ej — sj, and dropping the primes, the action becomes 



+ 2 7 i / dt [-Tiei(ei£i + 1/2) + 



j'+i _ £ j 



) - 2£ j £ 3 + l( £ j+l ~ £ j) ] 



+ 2 7L / dt [-T L (e L + s(L + l))((e L + s(L + l))e L + 1/2) + (e L + s(L + l))e L ] 



(63) 



which shows that the (e~ s< 3) = (e~ a ( L+1 >Q L ), where Ql is the time-integrated current flowing between site L and the right 
thermal bath. We shall denote by s' = (L + l)s. An additional change of fields, which leaves the bulk part invariant (see 0711 ). 
allows to further simplify the boundary terms: we set e' = x 2 \ Se and s' = \e(\ + 2ee), and we obtain (dropping the primes) 



J 3=1 J j=i 

+ 2 7 i J dt [-Tiei + eiei] 

- 7 l / dt [(s' + 2e L (l + s'e L )){T L + e L {s'T L - 1))] 



j+i - ej) - 2ejEj +1 (e j+1 - Ej) 



We differentiate S given in ( 164-b with respect to 



sr. 



VjE - 4Vj£ EjE j+1 



V,'_i£ — 4V 9 -_i£ Ea_iE, 



= 



where we used the notation WjX = Xj + \ — Xj. One thus has a constant K\ such that 



4 £ j £ i+i 

Differentiating now < f64b with respect to Ej one has 

VjE - Vj- 1 E+2(VjE) 2 Ej +1 + 2(V j -XE) 2 E j - 1 = 

and substituting d66l ) to get an equation on the Ej 's only, one obtains 

I- p.- -, . „ . 

e,-+ie,--i) = 



e 3+l £ j £ j-l 



The trick is to multiply d68t by 



which leads to 



£ j+i ~ £ j-i , 
£ i+i + £ j-i 



rx?-(v je ) a - 




r^-(v,_ l£ ) 2 


£ j+i £ i 




£ i £ i-i 



= 



and thus there exists a constant K 7 such that 



We substitute d66l l into the bulk part of the action j64l , and we arrive at 



-Stalk 



L - 1 



A" 2 



(64) 



(65) 



(66) 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 
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We differentiate the action with respect to the fields at the boundaries E\ , e\\ 

-(e 2 - £i) + 4(e 2 - £i)ei£2 + 27iei = 
~(e 2 - ei) - 2(e 2 - ei) 2 e 2 + 271 (ei -T x ) = 

Differentiating with respect to ex and el one gets 



(e L - £l-i) - 4(e L - £l_i)£l£L- 







2e L + (1 - 4T L e L + 4e L e L )A - T L (1 + 4e L e L )s' 2 
(ex - £ L -i) - 2(e L - £ L -i) 2 £l-i - 2 7L (1 + s'e L )(T L + e L {XT L - 1)) = 
We now proceed with solving the microscopic equations d66b and d7Tb - We search for a solution in the form 



tj = Asinh [2((j - l)B + C)] 



e j =E+—ta,nh[(j-l)B + C] 



(73) 
(74) 



(75) 
(76) 

(77) 



where A, B, C and E 1 are four constants to be determined by the four saddle point equations at the boundaries. We first note 
that, quite remarkably, d77b is an exact solution of the microscopic bulk saddle point equations d66b and ( P/ll . on condition that 



K\ = -A sinh(2B) 
One checks that the saddle equations at site 1 are solved by 

E = T, ; 



K 2 = sinh 2 B 



1 sinh 2B 

2 sinh 2C 



7i 



Eliminating 7^ between the saddle equations at site L yields 



A 2 = 



s' (T L s> - 1) 



16(T 1S > + l)(T L +T 1 (T L s> -I)) 



Substituting this result into ( T73T >. one gets 



sinh 



2((L - 1)B + C + e) 



— sinh 2C = 

1L 



where e is such that 

sinh 2 e = lo 

and oj is given by 

lo = \(T 1 -T l -XT 1 T l ) 
in accordance with ( fT2b . One can now eliminate C using j79l ) and this gives an equation involving B only 



sinh" 



2(L- l)B + 2e 



7r 2 + ll 2 + 2(7i7l) _1 cosh 2(L - l)B + 2e 



sinli 2B 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



(84) 



The large deviation function is given by the value of Wq l at saddle. Combining the bulk contribution d72| ) together with the 
boundary terms read from d64l i. one obtains 



1 

71 - 
L — 1 



1 



^(A) = | 7 i - -V7i 2 + sinh 2 2B + "- lL - \\f^+ sinh 2 2B 



■ sinlr B 



(85) 
z 

where B is solution of (l84l . 

Though the expressions are cumbersome, one can still solve the equation ( f84T > perturbatively in powers of L, by writing B in 
the form B = Bq/L + Bi/L 2 + . . .. To lowest order, B is of order 1/L and e of order 1: in ( |84l ), the right-hand term is negligible 
and one obtains Bq = — e = — arcsinh which yields the macroscopic fluctuation theory result /i(A) = — |(arcsinh t/uj) 2 
found previously in ( I9I12I I, as expected. To the next order, one gets 



V'qJA) = -/x(A) + 



1 - 5- — -J— ) M(A) 
271 27L 



+ 0(L- S ) 



(86) 
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which matches the announced result (f22t . 
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